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Abstract 

Exact expressions of the Pauli exclusion operator Q in the nuclear matter calcula- 
tion are presented in detail. Exact formulae are also given for the calculations of the 
single-particle-potential energy and the binding energy per nucleon with the exact 
Q operator. Numerical calculations of the G matrix in the lowest-order Brueckner 
theory are carried out to check the reliability of the standard angle-average approx- 
imation for the Q operator by employing the Bonn B and C NN potentials. It is 
observed that the exact treatment of the operator Q brings about non-negligible 
and attractive contributions to the binding energy. 
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1 Introduction 



It is one of the fundamental and open problems in nuclear structure physics 
to understand the saturation property of nuclear matter starting from a nucleon- 
nucleon {NN) interaction [1-3]. Because of the presence of strongly repulsive 
components in the short-range part of the A^A^ interaction, the nucleon- nucleon 
scattering correlation has a predominant importance in the nuclear matter 
calculation. The Bethe-Goldstone equation has been used for solving a two- 
nucleon scattering problem in nuclear matter. The influence of nuclear medium 
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on the nucleon-nucleon scattering is taken into account by considering the 
Pauh exchision principle and the self-energy effect on scattering micleons. 
Therefore the accurate treatment of the Pauli exclusion operator has been 
one of the essential requirements for the numerical calculation of nuclear mat- 
ter. 

The Pauli exclusion operator depends, in principle, not only on the magni- 
tudes of the total and relative momenta of scattering two nucleons but also 
on their angles. This angular dependence leads to couplings among partial 
waves, which makes numerical computations difficult. The difficulty due to 
the angular dependence has been avoided by employing the angle-average 
approximation[4,5]. The angle-average approximation, however, does not nec- 
essarily have quantitative justification, although there have been some studies 
to assess the reliability of this approximation[7-12]. 

One of the purposes of this work is to give analytic formulae for the cal- 
culation of the nuclear matter binding energy in the exact treatment of the 
Pauli exclusion operator Q. Furthermore we want to clarify how the G matrix 
depends on angular and linear momenta. With analytic expressions of the op- 
erator Q and the G matrix we perform numerical calculations to investigate 
the effect of the exact treatment of the operator Q on the nuclear matter 
binding energy by solving the Bethe-Goldstone equation. 

Very recently Schiller, Miither and Czerski[13] reported their calculations 
of nuclear matter properties using the exact Pauli operator Q. They pointed 
out that the standard angle-average approximation for the operator Q tends 
to underestimate the binding energy per nucleon at low densities but over- 
estimate it at higher densities. The exact treatment of the operator Q thus 
leads to a non-negligible improvement in the calculated saturation points. It 
would be interesting and deserving to examine whether the similar result can 
be reproduced in an alternative approach to the nuclear matter calculation 
with the exact Q operator. 

The organization of this paper is as follows. In Section 2 we give an exact and 
analytic expression of the operator Q in angular-momentum-coupling states. 
Section 3 is devoted to derive formulae which are necessary for the calculation 
of the single-particle-potential energy and the binding energy per nucleon in 
the exact treatment of the operator Q. We show, in Section 4, numerical 
results for the calculations with the exact and angle- averaged Q operators 
by employing the Bonn B and C A^A^ potentials[14]. We discuss the effect of 
the exact treatment of the operator Q and the accuracy of the angle-average 
approximation at various nuclear matter densities. Conclusions obtained in 
this study are given in Section 5. 
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2 Exact expression of the Pauli exclusion operator 



The reaction matrix {G matrix) in nuclear matter is defined by the Bethe- 
Goldstone equation 

G^v + v-G, (2.1) 
e 



where v is the two-body A^A^ interaction, e gives the starting energy minus 
the energy of an intermediate two-particle state, and Q stands for the Pauli 
exclusion operator. The operator Q prevents two particles from scattering into 
intermediate states with momenta below the Fermi momentum k-p, which is 
written as 

Q = ^ E Wl3){al3\Q{K - k^mh - ^f), (2.2) 

^ a/3 



where a is a single-particle state with the momentum k^-, the spin projec- 
tion Sa and the isospin projection Tq. The state \a(5) is a normalized and 
antisymmetrized two-nucleon state, and Q{x) the Heaviside step function. In 
the relative and center-of-mass (cm.) momentum system the operator Q is 
written as 

I K^dKdk I k^dkdk 

x\KkSaTaSpTp){{KkSaTaSf3Tf3\ - {K - kSpTpSaTa]} 

X e(|y + fe| - A;f) e(|y -k\- kp), (2.3) 

where lif and k are the cm. and relative angular momenta given by K—ka+kp 
and k—{ka—ki3)/2, respectively. 

We consider a partial wave decomposition of the relative state coupled to 
the angular momentum J with the cm. momentum K as 

\Kk{lS)JMTT,) =fisTj2i^^SS,\JM) 

m,Sz 

X I dkYim{k)\Kk)\SS,TT,), (2.4) 

where S is the total spin and I the orbital angular momentum of two-body rel- 
ative motion. We here have introduced the anti-symmetrization factor defined 
by fisT = [1 ~ (— 1)'^'^^"^]/2. The matrix element of the operator Q between 
angular-momentum-coupling states becomes 
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{Kk{hSi)JiM,T,T,,\Q\K'k'{l2S2)J2M2T2T,2) 

= 5{K - ^-SsiSi^TiTi^T.iT.i 

X g(/i JiMi, I2J2M2 : SiTikK9K<i>K), (2.5) 

where 



Q{kJiM^MJ2M2 : STkK0K(f)K) 



mim2Szi,Sz2 ' 
X {hmiSS.i I JiMi) {1217123 S,2\J2M2) 

X e(|y + fe| - M e(|y -k\- kp). (2.6) 

Here and (px arc the polar angles of K. Due to the presence of the step 
function, the solid angle k is restricted in the integration, and therefore / and 
J are, in general, not conserved in the matrix element of the operator Q, while 
K, k, S, T and are conserved. 

If we employ a reference frame, referred to as the K system hereafter, in 
which K points in the direction of the z axis, the angle between K and k 
coincides with the colatitude ^ of fe. In this case the hmit on the integration 
with respect to the angle 9 in Eq.(2.6) is given by —xq < cosO < xq, where 



for A; < ^k^-K'^/A, 



£!/l±r^ for 

1 otherwise. 



for ^kl-KyA<k<kF + K/2, (2.7) 



Since the integration with respect to the longitude of the relative momentum 
k can be carried out in the K system, the matrix element of the exclusion 
operator, denoted by Qq, is written as 

Qo{hJiMi,l2J2M2: STkK) 

= fhST hsT 5m,M2 Y.(-^f^^"l^^l-^jMhm\m{Ji -M1J2M1ILO) 

L 

xW{hJil2J2;SL) J PL{x)dx, (2.8) 



where Pl{x) is the Legendre polynomial, and we have used the notation / = 
\/2l + 1. It is clear from the above expression of the operator Q that, in the 
K system, the projection M is conserved though the magnitude of J is not. It 
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might be useful to show that the matrix element Qq{IiJiM,12J2M : STkK) 
satisfies the closure relation 

777^ E QoikJiM, kJ2M : STkK)5j,j, = fi.sT XoSj.jA.h- (2-9) 

We next consider the matrix element of the operator Q in an arbitrary 
reference frame in which the direction of K does not coincide with the z 
axis. By making rotation of the reference frame through the Euler angles 

— 4'K-iP = Ok and 7 = 0, we obtain a general expression of the matrix 
element of Q as 

Q{hhM^MJ2M2 : STkKeK(t>K) 

= E D^m'{Ok: (t>K. 0)Di^:M'iOK: 4>K. 0)Qo{hJiM' , I2J2M' : STkK) 

M' 

L 

XO 

xYlm{9k, 4>K)WihJil2J2; SL) J PLix)dx, (2.10) 



where the function Dj^j^i^a, (3,j) is the Wigner L'-function[15] of the Euler 
angles (a, (3, 7). In the derivation of Eq.(2.10) we have used the fact that only 
even integers L are allowed due to the parity conservation. We remark here 
that the matrix element of the exclusion operator is factorized as 

Q{liJiMi,l2J2M2 : STkK0K(f)K) 
= QihJiMi, I2J2M2 : STkK9K<t>K = 0) x e^(A^2-Mi)0K (2.II) 

The above expression shows clearly how the matrix element of the operator 
Q depends on the angle 

By using the recurrence formulae for the Legendre polynomials we obtain 
an integral formula for an even integer L as 

XO ^ 

j PL{x)dx = ^^[P^+,(a;o) - Pl-i{xo)]. (2.12) 


With use of the above relation we may rewrite Eq.(2.10) as 

Q{hJ,M^,l2J2M2 : STkKdK(t)K) 

= fhST fl2ST{xQ5i^i^5j^j^5MxM2 
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L>0,L=even 

xyLM(^i^,0i^)W^aiJi^2J2;-5L)[P^+i(xo) - Pl-i(2:o)]}- (2.13) 



It might be useful to give another expression of Eq.(2.10) written as 



Q{hJiMi,kJ2M2 : STkKeK(t>K) 

— fhST fhST [xoSi^h^J-i^J^^MiM'i 
L>0,L=even 

xYlm^Ok, <Pk)W{IMJ2; SL) ^^^^ {xl - l)Pl(a;o)], (2.14) 

where Pi{xo) is the first derivative of Pl{x) at a; = xq. It is easy to see from 
the above expression of the operator Q that if = 1, the second term on the 
right-hand side of Eq.(2.14) vanishes. Therefore ii k > kp + K/2, in this case 
Xo = 1 as given in Eq.(2.7), the matrix element of the operator Q conserves 
/, J and M and becomes unity. We emphasize that Eqs. (2.13) and (2.14) do 
not include integrals or step functions anymore, and are analytic and general 
expressions of the matrix element of Q which are valid for any vector K. 

We discuss here the angle-average approximation for the operator Q, which 
was first introduced by Brueckner and Gammel[4] and has been adopted in 
many nuclear matter calculations. The operator Q in the angle-average ap- 
proximation is defined as an average over the angle of K as 



Q{liJiM^MJ2M2 : STkK) 

= ^ [dk QihJiMi, hJ2M2 : STkK9K(f>K). (2.15) 

In this integration the second term on the right-hand side of Eq.(2.13) or 
(2.14) vanishes, and we have 

Q{liJiMi, I2J2M2 : STkK) = fi^sT Xo5i^i^5j^j^5miM2- (2-16) 



It is obvious that the first term on the right-hand side of Eq.(2.13) or (2.14) 
gives the angle-averaged Q operator and the second term provides the correc- 
tion which comes from non-spherical characters of the operator Q. 
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3 One-body potential and ground-state energy 



The G matrix equation in the lowest-order Brueckner theory is given in 
angular-momentum-coupling states by 



{kihSJiMiT\G{uj,K0K(l)K)\k2l2SJ2M2T) 

= ( h hJiST\v\k2l2J2ST)5j^j^5MiM2 

+ [k'dk{k,hJiST\v\kl[J[ST)5j,j>dM,Mi 

X Q{l[J[M[l'^4M!^ : STkROKCpK) x ^ 



u - T{K, k) 

X {kl2S4M^T\G{u;,KeK(l)K)\k2l2SJ2M2T), (3.1) 
where lu is the starting energy and T{K, k) the kinetic energy defined as 

T{K,k) = — + - (3.2) 
Am m 



with the miclcon mass m. Wc adopt the conventional QTQ spectra for the 
intermediate energies in the present work. 

The G matrix is generally diagonal in the cm. momentum K, the spin 
S and the isospin T, but not in the angular momentum I and the coupled 
angular momentum J of relative states. Corresponding to the factorization of 
the matrix element of the operator Q as shown in Eq.(2.11), the G matrix is 
factorized as 



{kihS.hM^T\G{u,KeK(t)K)\k2l2SJ2M2T) 

= {kihSJiMiT\G{uj,KeK(/)K = 0)\k2l2SJ2M2T) x e^(^2-Mi)0K_ (3 3) 

This relation expresses explicitly how the G matrix depends on the angle (px- 
It is to be noted that the G matrix is not invariant under the rotation about 
the z axis, but the rotation yields a phase factor as given in Eq.(3.3). 

The ground-state energy per nucleon of nuclear matter is given by 



A A//<pf 



^1) + ^ / u{k,)kldk,, (3.4) 



5'2m 

■ 
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where t is the one-body kinetic energy, u{kx) the potential energy of a nucleon 
in the occupied state |A), and pp the upper most occupied single-particle level 
(the Fermi level). Without loss of generality one can choose a reference frame 
in which the z axis (the quantization axis) coincides with the direction of the 
momentum kx as kx — (0, 0,/ca). We call this frame the L system hereafter. 
The self-consistent potential u{kx) is written in terms of the G matrix in the 
L system as 



^^(^a)= E(V|G'|A/x) 

(2T + 1) f f f 

= 2 J ^l^^i^ J ^"^o^^M J 

hJih-hST 0-10 
■m-im2SzMiM2 

x{kx^{hS)J,M,T\G{uj,KeKM\kx^{l2S)J2M,T){hmiSS,\J,M,) 
x{l2m2SS,\J2M^)Yi^^,{ex,Ax,)yilm2{K^K)^ (3-5) 

where uj is the starting energy given by 

u = {\\t\\) + u{kx) + {ii\t\iJi) + u{k^). (3.6) 



Here the angles (6*^, (p^) and iOxfj.) 0A/i) are the polar angles of A;^ and kx^ 
in the L system, respectively, where kx^i is the relative momentum defined by 
kxp^ = ikx-k^)/2. 

The integration in Eq.(3.5) is much simplified if we use the G matrix in the K 
system in which = = 0. Let {kx^,{hS)JiMT\Go{uj,K)\kxf,{l2S)J2MT) 
be the matrix element of the G matrix in the K system. Rotating the reference 
frame through the Eulcr angles {a = (j)^, /? = 6k, 7 = 0), the G matrix in the 
L system is related to that in the K system as 

{kx^,(hS)JiM,T\G(u;,Kek(t>K)\kx^(l2S)J2M2T) 
= E ^m,m(0k, Ok, 0)D'Xm{<I>k, Ok, 0) 

M 

X {kx^{hS)J,MT\Go{uJ,K)\kx^{hS)J2MT). (3.7) 



Substituting Eq.(3.7) into Eq.(3.5) and making integration with respect to 
the angle 4>^ we have 



hJil2J2MST -1 

x{kx,{liS)J^MT\G^{u,K)\kx,{l2S)J2MT) 



xF{kxkx^K{hJ,){l2J2)STM), (3.8) 
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z 




X 

Fig. 1. The relation among relevant momentum vectors. 

where we have used the facts that (pK — (f>fj, and (px^ = 0^ + 7r in the K system. 
The coefficient F in Eq.(3.8) is given by 

F(kxkx^,K(hJ,)(l2J2)STM) 

= Y: i-ir^-'^^D'^^^iO, 0K, 0)D'^MO, 9k, 0) 

mim2SzMiM2 

y<{limrSSz\JiM,){hm2SSz\J2M2)Yi^rn^{ex^,^)Yil^^^^^^^ (3.9) 
Using some formulae of the Wigner D-function, the coefficient F becomes 

F{kxkx^K{hJ^){hJ2)STM) 

= ^ ^(-1)^+"^ iiUih W{hJil2J2;SL){Ji -MJ2M\LQ) 

x{h0l20\L0) Pl{cos{9k + ^am))- (3-10) 



As shown in Fig.l, 9k + 9\fj, is the angle between the relative momentum 
k\fjt and the cm. momentum which satisfies 

cos^^ + 9x^) = ■ 3.11 



It is clear from Eq.(3.8) that the single-particle potential u(k\) is determined 
by the G matrix in the K system and the F coefficient which is a function of 
the angle 9k + ^A// between fc and K. 
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In the angle- average approximation the G matrix in the K system is inde- 
pendent of the quantum number M. Using the following closure relation of 
the coefficient F 



Y,F{kxkx^K{hJ,){kJ2)STM)6j,j, 



M 



2Ji + l 
An 




(3.12) 



the single-particle potential u{kx) with the angle-averaged Q operator becomes 



This expression of u{kx) agrees with the usual formula in the angle-average 
approximation for the operator Q. 



4 Numerical calculation 

In order to examine the effect of the exact treatment of the Pauli exclusion 
operator we performed a numerical calculation of the ground-state properties 
of nuclear matter by adopting the Bonn B and C A/^A/^potentials[14]. We solved 
self-consistently the coupled equation (3.1) for the G matrix in the K system 
where the cm. momentum K points in the z direction. In the calculation 
we took into consideration rigorously the contributions of the partial waves 
of J < 6. Other higher partial waves were taken up to J = 18 in the Born 
approximation. We checked the stability of the calculated result with respect 
to the number of mesh points in numerical integration. Furthermore we con- 
firmed that the same result was obtained within numerical errors by using two 
computer codes made independently. 

The calculated results of the binding energies per nucleon are presented in 
Table 1 and shown in Figs. 2 and 3, respectively. From Table 1 it is seen, 
as a common characteristic of the results for two NN potentials, that the 
exact treatment of the operator Q brings about attractive contributions to 
the binding energy per nucleon at any nuclear densities, compared with the 
result in the standard angle-average approximation. As is shown in Figs. 2 
and 3 the saturation densities hardly change in the exact treatment of the 
operator Q. 




x{kx^{lS)JT\Go{u;,K)\kx^{lS)JT). 



(3.13) 
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Table 1. Calculated results of the binding energy per nuclcon with the exact and 
angle-averaged Q operators for the Bonn B and C AW potentials[14]. The energies 
are in MeV. 



(1) Bonn B NN potential. 



A:F[fm-i] 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


exact 


-10.28 


-11.72 


-12.98 


-13.95 


-14.47 


-14.42 


-13.61 


average 


-10.18 


-11.59 


-12.82 


-13.75 


-14.25 


-14.16 


-13.32 


difference 


-0.10 


-0.13 


-0.16 


-0.20 


-0.22 


-0.26 


-0.29 


(2) Bonn C NN potential. 


fcF[fm-^] 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


exact 


-9.62 


-10.82 


-11.78 


-12.37 


-12.46 


-11.90 


-10.53 


average 


-9.52 


-10.69 


-11.62 


-12.18 


-12.24 


-11.65 


-10.24 


difference 


-0.10 


-0.13 


-0.16 


-0.19 


-0.22 


-0.25 


-0.29 



E/A [MeV] 




1.2 1.4 1.6 1.8 
kp [fm-^] 



Fig. 2. Calculated binding energies 
per nucleon as a function of the 
Fermi momentum for the Bonn 
B potential. 



E/A [MeV] 




1.2 1.4 1.6 1.8 

kp [fm-^] 



Fig. 3. Calculated binding energies 
per nucleon as a function of the 
Fermi momentum for the Bonn 
C potential. 
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These observations in the present work are quite different from those ob- 
tained in a very recent study by Schiller, Miither and Czerski[13]. Their results 
show that the exact treatment of the operator Q brings about repulsive contri- 
butions at higher densities, while small attractive contributions were obtained 
at low densities. This characteristic density dependence of the effect from the 
exact treatment of the operator Q leads to a non-negligible shift in the calcu- 
lated saturation point. 

In order to analyse the features of our calculated results in more detail, 
we introduce an approximate decomposition of partial waves in the exact 
treatment of the operator Q. Since the G matrix in Eq.(3.8) in the K system, 
Gq{uj, K), is not diagonal in / and J, as was discussed in the preceding section, 
the partial wave decomposition is not possible in a usual sense in the exact 
treatment of the operator Q. However, even in this case, we may decompose the 
contributions to the total potential energy per nucleon into two groups of terms 
which are diagonal and non-diagonal in / and J. We expect that the partial 
wave decomposition can be physically meaningful when the magnitudes of the 
contributions of the / J- non diagonal terms are sufficiently small in comparison 
with the ones of the Z J-diagonal terms. We call the diagonal term in I and J 
the partial wave contribution of the channel {S, T,l, J). The summation over 
M is taken in each of the partial wave contributions. 

We show in Table 2 the partial wave contributions of the channels with lower 
I and J, the sum of the Z J-nondiagonal terms and the higher partial wave 
contributions which are treated in the Born approximation. The contributions 
of the / J-nondiagonal terms are shown to be indeed small in the numerical 
calculation. The net contributions of the terms non-diagonal in / and J are 
estimated to be at most 0.4% of the total potential energy. Therefore we 
may say that the partial wave analysis could be meaningful even in the exact 
treatment of the operator Q. 

We observed from Table 2 that about sixty percent of the difference between 
the exact and approximate treatments of the operator Q comes from the IJ- 
nondiagonal terms and the remainder does from the Z J-diagonal terms. We also 
note that the P waves give the most important and attractive contributions. 
The S waves give the secondary important and repulsive contributions to the 
difference. The difference becomes, in general, smaller as / (or J) increases. 
The reason is that the Born approximation works well in higher partial waves 
and the G matrix approaches to the bare NN potential, which is independent 
of the operator Q. 
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Table 2. Partial wave contributions to the binding energy per nucleon with the 
exact and angle-averaged Q operators for the Bonn B iVA'^ potential[14]. The row 
of "non-diag" means the sum of the i J-nondiagonal terms. The row of "higher" is 
the sum of contributions of higher partial waves of 7 < J < 18 treated in the Born 
approximation. The energies are in MeV. 



1 rr — 1 ^ 

A;F[fm ] 


1.2 


1.2 


1.5 


1.5 


1.8 


1.8 


channel 


exact 


average 


exact 


average 


exact 


average 


1 /-f 


-13.018 


-13.019 


-20.112 


-20.117 


-26.749 


-26.760 




-1.419 


-1.419 


-3.673 


-3.671 


-7.741 


-7.737 




-0.227 


-0.227 


-0.708 


-0.708 


-1.674 


-1.674 


'h 


-0.049 


-0.049 


-0.195 


-0.195 


-0.540 


-0.540 




-15.638 


-15.637 


—22.195 


—22.210 


-28.058 


-28.109 




0.909 


0.910 


2.300 


2.301 


4.605 


4.607 




-2.393 


-2.393 


-5.933 


-5.930 


-11.813 


-11.808 




0.145 


0.147 


0.426 


0.430 


0.920 


0.932 




0.105 


0.105 


0.365 


0.365 


0.931 


0.931 




-0.369 


-0.369 


-1.212 


-1.212 


-2.957 


-2.957 


^G5 


0.038 


0.038 


0.150 


0.150 


0.418 


0.419 




0.016 


0.016 


0.072 


0.072 


0.222 


0.222 




-0.071 


-0.071 


-0.301 


-0.300 


-0.873 


-0.873 




2.982 


2.989 


6.443 


6.462 


11.944 


11.987 




0.492 


0.493 


1.320 


1.320 


2.783 


2.784 


1 TT 


0.104 


0.104 


0.363 


0.363 


0.911 


0.911 




-2.552 


-2.552 


-4.535 


-4.536 


-6.216 


-6.218 


3 7-) 


6.690 


6.693 


14.980 


14.988 


28.680 


28.702 




-4.423 


-4.415 


-10.859 


-10.839 


-21.637 


-21.598 




-0.305 


-0.305 


-0.908 


-0.908 


-2.015 


-2.014 


SF3 


0.876 


0.876 


2.486 


2.486 


5.458 


5.458 




-0.148 


-0.148 


-0.579 


-0.578 


-1.661 


-1.660 


'H, 


-0.040 


-0.040 


-0.160 


-0.160 


-0.457 


-0.457 


'H, 


0.159 


0.159 


0.590 


0.590 


1.564 


1.564 




-0.016 


-0.016 


-0.077 


-0.077 


-0.250 


-0.250 




-0.007 


-0.007 


-0.034 


-0.034 


-0.115 


-0.115 


higher 


0.044 


0.044 


0.204 


0.204 


0.624 


0.624 


non-diag 


-0.083 


0.0 


-0.156 


0.0 


-0.226 


0.0 


Total 


-28.197 


-28.094 


-41.940 


-41.747 


-53.921 


-53.630 
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Table 3. Single-particle potential energies for various kx values at fcp = 1.40[fm~ ] 
in the exact and angle-average treatments of the operator Q for the Bonn B NN 
potential [14]. The other notations are the same as in Table 1. 



kx[im-'] 


u{kx) 
exact 


average 


difference 


0.047 


-87.96 


-88.54 


0.58 


0.237 


-87.33 


-87.84 


0.51 


0.533 


-84.69 


-84.90 


0.21 


0.867 


-79.33 


-79.12 


-0.21 


1.163 


-72.69 


-72.22 


-0.47 


1.353 


-67.71 


-67.22 


-0.49 



Next, we show in Table 3 the single-particle potentials in the exact and 
angle-average treatments of the operator Q and their differences for the val- 
ues of kx of the mesh points in the Gaussian integration. The differences of 
the single-particle-potential energies in the exact treatment from those in the 
angle-average approximation are found to be rather small for all the k\ values. 
It is remarkable to see that the angle-average approximation leads to overes- 
timation of the single-particle-potential energy for smaller values of kx and 

A u [MeV] 



T 1 1 1 1 1 1 1 r 




0.5 1 



Fig. 4. Differences of n(A;;^)'s in the exact and angle-average treatments of the op- 
erator Q as a function of kx/kp at various nuclear matter densities for the Bonn B 
potential. The quantity Ait is defined as the difference of u{kx) with the exact Q 
from that with the angle-averaged Q. 
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does to underestimation for larger values of k\. This feature is observed at 
any nuclear matter densities both for the Bonn B and C potentials. We dis- 
play in Fig. 4 the kx dependence of the differences between two m(A;a)'s in the 
exact and angle-average treatments of the operator Q at three nuclear matter 
densities. As shown in Table 3 and Fig. 4 the differences change the sign at 
about k\ = kp/^. This fact implies that the angle-average approximation gives 
the correct result at around k\ = kp/2. Considering larger phase volumes for 
larger values of kx, it is understandable that the angle-average approximation 
tends to underestimate the binding energy per nucleon. 

We may conclude from the present study that the non-spherical character 
of the operator Q causes a non-negligible and attractive effect on the nuclear 
matter binding energy per nucleon, which has been disregarded in the standard 
angle-average approximation for Q. We believe that the analytic expression 
given explicitly in this paper for the operator Q in the G matrix is useful in 
making a precise numerical calculation of nuclear matter. 



5 Concluding remcirks 

We derived an exact and analytic expression of the matrix element of the 
Pauli exclusion operator in nuclear matter, and gave its relation to that in 
the angle-average approximation. Furthermore we discussed the breaking of 
relevant symmetries in the matrix element of the exclusion operator. We pre- 
sented the rigorous expressions of the single-particle-potential energies and the 
ground-state energy in a form that would be suitable for a precise numerical 
calculation in the exact treatment of the operator Q. 

In order to examine the effect of the exact treatment of the operator Q and 
assess the reliability of the angle- average approximation for the operator Q 
we performed numerical calculations of the ground-state properties of nuclear 
matter for various kp values by employing the Bonn B and C NN potentials. 
We found that the exact treatment of the operator Q brought about non- 
negligible and attractive contributions to the nuclear matter binding energy 
per nucleon. Our calculations also clarified the degree to which the angle- 
average approximation is reliable quantitatively. 
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